Astrometric Microlensing of Distant Sources due to the Stars in 

the Galaxy 



Mareki Honma 1,2 & Tomoharu Kurayama 1 ' 3 



ABSTRACT 

We investigate properties of astrometric microlensing of distant sources (such 
as QSOs and radio galaxies) caused by stars in the Galaxy, mainly focusing 
on application to the VERA (VLBI Exploration of Radio Astrometry) project. 
Assuming typical parameters for the Galaxy disk and bulge, we show that the 
maximum optical depth for astrometric shift of 10 //as-level is 8.9 x 10~ 2 for 
QSO-disk lensing case and 3.8 x 10~ 2 for QSO-bulge lensing case. We also find 
that the maximum optical depth for QSO-disk lensing is larger by an order 
of magnitude than that for disk-disk or bulge-disk lensing case (assuming a 
typical source distance of 8 to 10 kpc). In addition to optical depth, we also 
calculate the event rate and find that the maximum event rate for QSO-disk 
lensing case is 1.2xl0 -2 event per year, which is about 30 times greater than 
that for disk-disk lensing case. This high event rate implies that if one monitors 
10 QSOs behind the Galactic center region for 10 years, at least one astrometric 
microlensing event should be detected. Therefore, monitoring distant radio 
sources with VERA can be a new tool to study astrometric microlensing caused 
by stars in the Galaxy. 

We also study the event duration of astrometric microlensing, and find that 
the mean event duration for QSO-disk lensing is 7.5 yr for QSOs located near 
the galactic center. This event duration for QSO-disk lensing is reasonably short 
compared to the project lifetime of VERA, which is anticipated to be ~ 20 yr. 
We also find that while the minimum event duration for bulge-bulge lensing is 
as short as 2.6 yr, the event duration for disk-disk lensing case cannot be shorter 
than 15 yr. Thus, although astrometric microlensing of bulge sources/lenses can 
be studied by optical astrometric missions like SIMA and GAIA, detections of 
disk events with the space astrometric missions are fairly difficult because of the 
limited project lifetime (typically ~5 yr) as well as the heavy dust extinction. 
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Therefore, for studying astrometric microlensing by disk stars, VERA can be 
one of powerful tools by observing distant sources like QSOs and radio galaxies. 
We discuss the implications of astrometric microlensing for VERA by focusing 
on estimating the lens mass, and also present some possible candidates of radio 
sources toward which astrometric microlensing events should be searched for 
with VERA. 

Subject headings: Galaxy : disk — Galaxy : structure — gravitational lensing 
— stars : low-mass, brown dwarfs 



Gravitational microlensing is one of the most promising tools to study invisible lenses 
like MACHOs and low-mass stars in the Galaxy. Paczynski (1986) first proposed to use 
microlensing effect to search for MACHOs in the Galaxy's halo based on a photometric 
monitoring of millions of stars in the Galaxy's bulge and the Magellanic Clouds. A number 
of groups have been conducting such 'photometric' searches of microlensing to detect source 
magnification by gravitational microlensing (e.g., MACHO, EROS, OGLE, MOA, etc.), 
finding more than a hundred of photometric microlensing events. 

In addition to such photometric microlensing, there is another type of microlensing 
called 'astrometric microlensing', in which the positional shift of the lensed image, rather 
than the magnification of the source, is used to detect microlensing events. Recently, a 
number of studies have been made on such 'astrometric' microlensing events (e.g., Hosokawa 
et al. 1993; 1997; H0g, Novikov k Polnarev 1995; Miyamoto k Yoshii 1995; Walker 1995; 
Paczynski 1996; 1998; Miralda-Escude 1996; Boden et al. 1998; Dominik k Sahu 2000). 
One of the major findings of these studies is that the probability of astrometric microlensing 
is much larger than that of photometric microlensing (e.g., Miralda-Escude 1996; Hosokawa 
et al.1997; Dominik k Sahu 2000; Honma 2001). For instance, the size of a photometric 
lens is equal to the Einstein-ring size given by 



Here M is the lens mass, and D d , D ds , and D s are the observer-lens distance, lens-source 
distance, and observer-source distance, respectively. If a source comes within this radius 
from the lens, a source is magnified by more than a factor of 1.34. Meanwhile, the size of an 
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astrometric lens is larger than R E by a factor of /9 max , which is given by following equation 
(e.g., Honma 2001). 

Anax = #E/#min- (2) 

Here #e is the angular size of the Einstein-ring radius (6>e = Re/D^), and 6> min is the 
minimum angular shift that can be detected by astrometric observation. Thus, if a source 
comes within /3 max i?E from the lens, the source position is shifted by larger than # m i n . For 
astrometric lensing, /5 max can be as large as 50 ~ 100 assuming 9 min of 10-/xas level (e.g., 
Honma 2001). 

Such a high astrometric accuracy, although not yet achieved, will be available soon 
because a number of astrometric missions are planned in early 21st centuries. For instance, 
there are four space astrometric missions to be launched by ~2010: DIVA, FAME, SIM 
and GAIA. Astrometric accuracies anticipated for those missions are: 200 /ias for DIVA, 50 
/xas for FAME, and 10 /zas (or possibly higher) for SIM and GAIA. In addition to those 
space missions, there is another ground-based astrometric mission called VERA (VLBI 
Exploration of Radio Astrometry, e.g., Sasao 1996; Honma, Kawaguchi & Sasao 2000 and 
reference therein), which utilize a phase- referencing VLBI technique for astrometry of radio 
sources. The most remarkable difference between VERA and space astrometric missions is 
that VERA can observe thousands of distant sources like QSOs and radio galaxies to trace 
effect of astrometric microlensing. An advantages of using distant radio sources is that the 
column density of lens can be much higher, leading to higher event probability. 

Investigating disk stars by astrometric microlensing with missions described above has 
significant scientific merits. First, astrometric microlensing provide us a new tool to study 
the lower end of stellar mass function, which is thought to be dominant populations in the 
Galaxy's disk. To date, the shape of the mass function remains uncertain in particular at 
its lower end, because the lower main-sequence stars are too faint for optical observations. 
Studying the low-mass star may also inform us about the nature of dark matter in the 
disk, as the disk dark matter may be baryonic matter in the form of low-mass stars and/or 
brown dwarfs. Further, studies of astrometric microlensing by disk stars will also provide 
fundamental information on parameters of the Galaxy's disk. For instance, one can estimate 
the disk scale length from the optical depth distribution with the Galactic latitude / (e.g., 
Dominik & Sahu 2000). Also, it is possible to extract the disk density from astrometric 
microlensing events provided sufficient number of events are observed for statistics. From 
the scale length and disk density, one can estimate the total mass of the Galaxy's disk, and 
thus the density profile of dark halo inside in the disk region can be also deduced. 

While the astrometric microlensing of stars are extensively studied in previous studies 
focusing on implications for SIM and GAIA(e.g., Dominik & Sahu 2000), there has been 
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no study on astrometric microlensing of distant radio sources by disk stars considering 
an application to VERA. The most remarkable aspect of using distant sources is that 
some important parameters such as the lens mass may be determined independently of 
lens distance (e.g., Honma 2001), as well as high event probability due to large column 
density toward a source. For these reasons, in the present paper we extensively study the 
astrometric microlensing of distant sources by disk stars, and discuss its implications for 
VERA. 

The plan of this paper is as follows: in section 2 we present calculations of optical 
depth, and in section 3 event rate. We also compare results for QSO-disk/bulge lensing 
cases with disk/bulge-disk lensing cases. In section 4, we present event durations of 
astrometric events, and in section 5 we discuss how the lens mass can be estimated from 
event durations. Finally in section 6, we discuss the implication of astrometric microlensing 
for VERA, and also summarize the findings of the present paper. 



2. Optical Depth 

2.1. basic equations 

The optical depth of astrometric microlensing event is defined as follows (e.g., 
Miralda-Escude 1996): 

r as t = / dD d . (3) 

Here, Re is the Einstein-ring radius, p is the mass density in the lens, is the lens 
distance from the observer, and M is the lens mass. Note that this lens mass should be 
regarded as the mean lens mass averaged over the possible range of stellar mass function 
(the effect of stellar mass function will be discussed later in the present paper). The factor 
Anax, defined in equation (|), denotes the ratio of astrometric lens size to the Einstein ring 
size. The factor of /^ ax in equation (^) is crucial to the difference of the optical depth 
between astrometric and photometric microlensing. As is mentioned earlier, in the distant 
source case /? max could be as large as 50 ~ 100, and thus the optical depth for astrometric 
microlensing could be larger than that of photometric microlensing by 3 to 4 orders of 
magnitudes. By substituting the expression of Re and /9 max into equation @, one can 
obtain following expression of the optical depth, 

T - = ^~/ p a) dDi - (4) 

Interestingly, the optical depth is proportional to the mean lens mass M, in contrast to 
photometric microlensing in which the optical depth is independent of the mean lens mass 
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M. 



2.2. mass model 

For the density distribution in the Galaxy's disk, we assume an exponential disk for 
both radial and vertical profiles, namely, 

p d (R, z) = po exp d R ° - . (5) 

Here, po is the disk density in the vicinity of the Sun, and d and h are the radial and vertical 
scale lengths, respectively. We take rather conservative values for these parameters as 
Po = 0.08M Q pc~ 3 , d = 3.5 kpc, and h = 300 pc (table 1 summarizes the model parameters 
used in this paper). Note that the disk model and its parameters considered here are the 
same as those used in previous studies (e.g., Hosokawa et al. 1997; Dominik, Sahu 2000; 
Honma 2001). 

In order to take into account the contributions from bulge stars, we consider a bulge 
model whose density distribution is given by a spherical Plummer model. The density 
distribution of the bulge is given by 

where M b is the total mass and a is the scale length of the bulge, respectively. Although 
the bulge may be elongated and bar-like system as suggested by recent studies, here we do 
not consider such an asymmetric density distribution, as the precise parameters of bulge 
elongation is still highly uncertain. For the bulge parameters, we assume the bulge mass 
M b of 0.8 x 10 10 M Q and the scale length a of 1 kpc (e.g., Dehnen & Binney 1998). 

Here we also introduce the stellar mass function to calculate the mean lens mass M. 
To be conservative, we assume the Salpeter-type stellar mass function with the power law 
index of 2.35, namely, 

4>(m) = C(f,m~ 2 ' 35 (for my < m < my). (7) 

Here gives a scaling constant, and and my are the lower and upper cutoff of the 
stellar mass function. Based on that mass function, we obtain the mean lens mass as 

M = I m<P(m) dm 
J (j)(m) dm 
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With a lower cutoff M L of O.O8M , which corresponds to the hydrogen burning limit, 
we obtain the mean mass of 0.3M Q (note that M is independent of M\j provided My is 
sufficiently large). This mass M of O.3M corresponds to that of lower main-sequence stars, 
which are thought to be dominant populations in the Galaxy's stellar component. Thus, in 
the following sections we assume the mean lens mass of O.3M except otherwise noted. Of 
course, the mean mass is quite uncertain since we have little knowledge of the lower end of 
the stellar mass function. In case of different mean mass, the optical depth for mean lens 
mass M' can be obtained by multiplying a factor of (M'/M) to the optical depth obtained 
in the present paper. 

2.3. results 

Using equations described above, one can obtain the optical depth for astrometric 
microlensing. In the present paper, we consider several cases with different source and lens 
populations. Table 2 summarizes the cases we consider in the present paper. For QSO-disk 
and QSO-bulge lensing, in which distant QSO is being lensed by disk or bulge stars, the 
distance to QSO is assumed to be infinite. On the other hand, for disk-disk lensing case, 
we assume the source distance D s of 8 kpc. The first reason for setting D s = 8 kpc is that 
the distance corresponds to that to the disk stars near the Galactic center, where the disk 
density is highest. The second reason is that bright stars like AGB have My of mag, 
and hence their apparent magnitude becomes about 15 mag at the distance of 8 kpc in 
case of no interstellar extinction. Since the stars at low galactic latitude significantly suffer 
from strong absorption by the interstellar medium (later we will show that we have to 
observe stars at fairly low latitude to detect an astrometric microlensing), their apparent 
magnitude could be reduced to ~20 mag or even fainter, which is close to or below the 
limiting magnitude of optical astrometric missions like SIM and GAIA. Meanwhile, for 
bulge-disk or bulge-bulge lensing case, we assume a slightly longer distance, D s of 10 kpc. 
This is because a part of bulge sources can be seen through low-extinction regions like the 
Baade window, which are a few degrees away from the Galactic center. In these windows, 
we may be able to observe slightly far sources compared to disk sources partly because the 
extinction is low, and partly because there are quite many bright stars like bulge red giants 
(note that these 'windows' are not suitable for search for disk lens event since their galactic 
latitude is too high for disk lensing; for instance, b ~ 4° for the Baade window). 

Figure 1 shows the optical depths calculated for the cases described above. First of all, 
for QSO-disk/bulge lensing case, the maximum optical depth is obtained for sources behind 
the Galactic center (i.e., (Z, b) = (0°, 0°)), being 8.9 x 1(T 2 for disk lens, and 3.8 x 1(T 2 
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for bulge lens, respectively. This result indicates that nearly one of ten sources in this 
direction is always lensed by disk stars, and the optical depth of bulge lens is nearly half of 
that of disk lens (we will discuss how to discriminate disk lens event from bulge lens event 
later). However, the optical depth for bulge lens rapidly decrease with increasing I, and its 
contribution becomes negligible at I > 10°. For QSO-disk lensing, the dependence on b is 
much stronger that on that I. For instance, the optical depth at (I, b) = (0°, 5°) is smaller 
than that at (Z, b) = (0°, 0°) by almost an order of magnitude. This strong dependence 
comes from the short vertical scale length h (300 pc) compared to the radial scale length d 
(3.5 kpc). These results indicate that one should look for astrometric microlensing by disk 
stars at fairly low Galactic latitude (e.g., less than 2 degrees). 

For comparison, we also show in figure 1 the optical depth for disk-disk lensing case 
and bulge-bulge/disk lensing cases. From figure 1 one can find that the maximum optical 
depth for disk-disk lensing is 6.2xl0~ 3 , being smaller by a factor of 15 than the maximum 
optical depth for QSO-disk lensing. Also, the maximum optical depth for bulge-disk lensing 
case is 9.3xl0~ 3 , 10 times smaller than that for QSO-disk lensing. The optical depth for 
bulge-bulge lensing case is far smaller, being l.OxlO -3 . These significant differences of 
optical depths between distant source and disk/bulge source cases can be understood if 
one carefully looks at the equation (f|). In the distant source case, the factor (1 — D&/ D s ) 2 
is almost unity within the entire Galaxy, and thus the largest contribution to the optical 
depth occurs at the galactic center, where the disk density is highest. On the other hand, 
in the disk-disk lensing case the largest contribution to the optical depth comes from stars 
near the Sun (where the disk density is relatively low), because the factor of (1 — D&/ ' D s ) 2 
becomes fairly small for lenses near the Galactic center. This fact also explains that for the 
disk-disk lensing case the dependence of r ast on I is not so strong as that for the QSO-disk 
case: the density of local stars, which contribute most to the optical depth, does not vary 
strongly with changing I. 



3. Event Rate 

3.1. basic equation 

The event rate is another important quantity in microlensing study, since the event 
rate is directly related to the observational strategy. The event rate is usually defined 
as the number of lenses coming inside the microlensing tube per unit time (e.g., Griest 
1991). The radius of the microlensing tube is given by the lens size, which is /? max -RE for 
astrometric microlensing. By definition, the event rate strongly depends on the distribution 
of tangential velocity of lenses, and thus is usually a function of velocity distribution of lens 
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population. In this section, for simplicity, we assume that the tangential velocity of lens 
relative to the tube v± is uniquely determined at a given distance D d . This v± corresponds 
to an averaged tangential velocity at a given distance, which can be obtained by integrating 
the velocity distribution function over the possible velocity range (in later section we will 
treat the effect of velocity distribution function more precisely). Under this assumption, 
the event rate for astrometric microlensing T ast can be written as 

r as t = J 2(3 max R E -^v ± dD d . (9) 

This expression is correct regardless of the shape of velocity distribution function. For 
instance, in case that the lens population has an isotropic velocity distribution, the number 
flux of lenses which crosses unit area per unit time (F) is given by 

71 M 

Note that only the population coming in the tube is considered but the population going 
out of the tube is ignored. Since the area of the lensing tube per unit length dD d is 
27r/9 max i?E^-Dd, one can obtain the event rate per unit length (dT^/ dD d ) as 

= 27r/3 maxJ R E x-^-v ± = 2/5 maxJ R E -j- vj_, (11) 
dD d 7T M M 

which is a differential form of equation (||). On the other hand, if the lens population is 

perfectly streaming in a direction without any random motion, the number flux of lenses per 

unit area per time is given by p/Mv± (assuming that the streaming motion is perpendicular 

to the unit area considered). Since the geometric cross section of the tube per unit length 

seen from the lens population is 2f3 max R^dD d , one may obtain the same expression for the 

differential event rate dT^/dD^ for the isotropic case. 

Substituting equations ([!]) and (§) into @, one may obtain the following expression of 
the event rate, 

Tast = ^-lpvJ^] dD d . (12) 



c 2 6 ■ J \ D 

Note that the event rate for astrometric microlensing is independent of the lens mass M, as 
is pointed out in previous studies (e.g., Dominik & Sahu 2000; Honma 2001). 



3.2. velocity model 

The tangential velocity of the lens relative to the tube can be defined using the motion 
of source, lens and observer as follows, 

v± = -p- v o± + - v d ±. (13) 
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Here v Q ±, v s ±, and v&± are the two dimensional velocities of the observer, source, and 
lens perpendicular to the line of sight. In order to take into account the effect of velocity 
dispersion in a rather simple manner, here we calculate the lens tangential velocity v± as 



v± = \l \v ± \ 2 + (7, 2 ens + ( ^ ) <r s 2 ource , (14) 



where o"i cns and o" sourcc are the velocity dispersion of the lens and source population, 
respectively. 

The velocities of disk stars can be easily described based on the so-called flat-rotation 
model. Since the rotation velocity of the Galaxy's disk is known to be almost constant at 
any galacto-centric distance, here we assume a constant circular velocity t> C i rc of 220 km/s 
(e.g., Kerr & Lynden-Bell 1986). Given the circular velocity in the disk and star position (/, 
b and the distance to the star), one can determine its tangential motion with respect to the 
line of sight caused by the galactic rotation. In addition to circular rotation velocity, stars 
in the disk are known to have velocity dispersion of 15 ~ 20 km/s, and thus in case of disk 
stars (acting as both lens and source) we set <7i ens and/or o" sourcc of 20 km/s. 

In figure 2 we show the distribution of v± with the lens distance Figure 2 clearly 
demonstrates the difference of tangential velocities in two cases; QSO-disk lensing and 
disk-disk lensing case. For the QSO-disk lensing case, the largest tangential velocity is 
obtained toward a source behind the Galaxy's center, which is as large as 440 km/s. This 
velocity, corresponding to 2xt> circ , is obtained because behind the Galactic center the lens 
motion becomes completely opposite to the observer motion. On the other hand, for the 
disk source case the largest tangential velocity is not larger than 100 km/s. This is due to 
the fact that the source is not further than the Galaxy's center, and thus the source as well 
as the lens are moving in the same direction as that of the observer. In particular, in the 
vicinity of the Sun (i.e., within 1 kpc from the Sun), where the contribution to the optical 
depth becomes the maximum for disk-disk lensing case, the tangential velocity v± is no 
larger than 40 km/s, which is an order of magnitude smaller than the largest tangential 
velocity for QSO-disk lensing case. These facts imply that the event rate for disk source 
case could be significantly smaller than that for distant source case. 

On the other hand, for the bulge model we use the isothermal velocity dispersion 
obtained from the Plummer potential. For the tangential velocity dispersion, the Plummer 
model gives (e.g., Lynden-Bell 1962) 

^ = 1 GM h 
3 \Jr 2 + a 2 

The parameters M b and a are the bulge mass and the scale length, and taken to be 
M b = 0.8 x 10 10 M o and a = 1 kpc, as described in the previous section. With these 
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parameters, the tangential velocity dispersion at the Galactic center becomes ~ llOkm/s, 
which is in agreement with observed velocity dispersion (e.g., Tiede Sz Terndrup 1999). The 
bulge is also known to have a rigid rotation with an angular velocity of ~ 50 km/s/kpc 
(Tiede & Terndrup 1999). In the present paper, however, we do not include this rigid 
rotation, since at the central bulge (where most of bulge lensing events occur), the tangential 
velocity is dominated by the velocity dispersion rather than the rotational motion. 

3.3. results 

Based on the tangential velocity and density distribution described above, one can 
calculate the event rate of astrometric microlensing using equation @. Figures 3 shows the 
event rate distributions for cases considered in the previous sections (see table 2). From 
figure 3 one sees that for QSO-disk lensing case the maximum event rate is 1.2 x 10~ 2 
event/yr, which is reasonably high for practical observations. In fact, this indicates that an 
event can be detected if one observes 100 sources for a year, or 10 sources for 10 years. If 
50 sources are monitored for the whole mission lifetime of VERA (~ 20yr), one can expect 
to find ~10 astrometric microlensing events. For QSO-bulge lensing case, the maximum 
event rate is 5.0 x 10~ 3 , nearly half of that for QSO-disk lensing. Thus, for QSO behind the 
galactic center, the total event rate will be around 1.7xl0~ 2 event per year. 

The fact that QSO-disk lensing and QSO-bulge lensing rates are comparable raises 
a new problem, that is how to discriminate disk lens events from bulge lens events. To 
do this, one can use an additional observable of astrometric microlensing event, which is 
the direction of lens motion. If the whole source motion along the circular trajectory is 
monitored, then one can directly determine the direction of lens motion; the lens motions is 
perpendicular to the line connecting between the original source position and the position 
of maximum image shift. For QSO-disk lensing, one can expect that the direction of the 
lens motion is along the galactic plane. For instance, a typical disk lens behind the galactic 
center is moving at 400 km/s with respect to the lensing tube, while that star is likely 
to have a random motion of only 20 km/s or so (disk velocity dispersion). Therefore, the 
alignment of the lens direction and the galactic plane should be within a few degrees (~ 
20/400 radian). On the other hand, since for QSO-bulge lensing case the lens tangential 
velocity is dominated by the random motion of bulge stars, there should be no tight 
correlations between the direction of the lens motion and the galactic plane. 

For disk-disk lensing case the largest event rate is 4.0 x 10 -4 , being ~ 30 times smaller 
than that for the QSO-disk lensing case, and also 10 times smaller than that for QSO-bulge 
case. At the galactic center (I, b) = (0°, 0°) the event rate for disk-disk lensing is only 
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1.3 x 10 4 . This small event rate comes from the small tangential velocity for the disk source 
case as well as the fact that stars around the Galactic center, where the stellar density is 
highest, make little contributions to the disk-disk lensing (as discussed for the optical depth 
calculation). At the galactic center the bulge-disk lensing case and bulge-bulge case give 
higher event rate than that for the disk-disk lensing case, being 1.3 x 10~ 3 and 4.0 x 10~ 4 , 
respectively. This is due to higher velocity dispersion in the bulge than in the disk. Hence, 
if one looks for astrometric microlensing events with space astrometric missions like SIM 
and GAIA, it is better to observe bulge sources rather than disk stars. 

Note that for disk-disk lensing cases the maximum event rate is obtained at / ~ 10° 
rather than toward the Galactic center. This is because the event rate (equation||) is 
obtained by integrating a product of the density p and tangential velocity v±: while the 
maximum density along the line of sight is decreasing with Z, the maximum tangential 
velocity is peaked around / ~ 10°. On the other hand, toward I = 0° the rotation velocities 
of the source, lens and observer cancel out, and the tangential velocity reduces to the 
velocity dispersion a, making the event rate fairly small. 



4. Event Duration 

4.1. basic equations 

The event duration, which is one of the most important observables in individual 
microlensing events, is defined as period during which a lens is located within a microlensing 
tube. For astrometric microlensing, the event duration can be described as 

, 2/3 max i?E c- 

*ast = COSC), (16) 

v± 

where 2/3 max i? E corresponds to the lens diameter and the angle S defines the lens-source 
impact parameter as (3 max R^ x sin 5. The averaged event duration can be obtained by 
integrating t ast over possible range of 5 as (e.g., Honma 2001) 

Substituting equations flU) and (§) into ([17[) , one can obtain 

2nGM /D As 



t 



c 2 v±6 min V D, 



Note that the average event duration is defined at any lens distance D d , and varies with 
changing D d through the factor of (D ds /D s ) as well as v±, which also depends on D d . Thus, 
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in order to obtain expected event duration toward a source (hereafter denoted as (t as t)), 
one has to integrate the average event duration t ast along the line of sight from observer 
to source. We also note that for distant source case, the averaged event duration becomes 
independent of the lens distance as (Dds/D s ) approaches to unity (Honma 2001). Hence, 
for distant source case, the uncertainty in lens mass determination can be smaller than that 
for disk source case. 

While one can obtain the expected event duration (t as t) by integrating equation (|18|) 
over possible range of Dd, here we use another simple way, that is to use the optical depth 
and event rate obtained in the previous sections. These two quantities are related to the 
expected event duration through the following equation, 

(*ast> = (19) 
J- ast 

Based on this equation together with the optical depth and event rate calculated in previous 
sections, we calculate the event duration. 



4.2. results 

Figure 4 shows the distribution of (t as t) with various lines of sight, for both of 
disk/bulge source and distant source cases. For QSO-disk lensing, (t ast ) is found to be 7.5 
yr toward the Galactic center (in case of b = 0°). Around the Galactic center region (e.g., 
/ < 30°), the expected event duration changes only slightly with I, but increases rapidly 
for the outer Galaxy lenses (I = 90° to 180°), reaching at 91 yr at the anti-galactic center. 
Hence, in order to detect events within rather practical time-scale (i.e., less than 20 yr), 
one should observe sources around the Galactic center regions. 

For QSO-disk lensing toward the galactic center region, the lenses may be divided into 
two populations based on their location along the line of sight, namely the lenses in front of 
the galactic center and the lenses behind the galactic center. For these two populations, the 
relative proper motions of the lens are completely different: while the proper motion of lens 
in front of the galactic center is dominated by the velocity dispersion in the disk, the proper 
motion of lens behind the galactic center is dominated by the galactic rotation (see figure 
2). Due to this difference, these two populations also have different distribution of event 
rate and event durations. To show this clearly, in figure 5 and 6 we show the dependence 
of differential event rate and event duration on the lens distance. As seen in figure 5, most 
of events are likely to be caused by the stars behind the galactic center because of the 
large proper motion. For instance, toward / = 0° about 95% of events are caused by the 
stars behind the galactic center. Figure 6 shows that the events caused by stars behind the 
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galactic center have shorter event duration because of larger proper motion. In fact, for the 
case of / = 0°, the lenses behind the galactic center gives a typical duration of 4.1 yr, while 
the lenses in front of the galactic center gives a duration of 91 yr. In practice, however, the 
distinction between two populations is not as clear as in figure 5 or 6 because of the effect 
of the mass function or velocity dispersion. For instance, a less massive star in front of the 
galactic center could cause shorter event, which may be comparable to events caused by 
more massive star behind the galactic center (in the next section, we will discuss the effect 
of the mass function and the velocity dispersion more precisely). 

On the other hand, for disk-disk lensing case, the dependence of (t ast ) on I (figure 4) 
is not simple, having multiple local maximum and minimum values. The minimum event 
duration is obtained toward I ~ 10° where the tangential velocity of the lens is largest (see 
figure 2). The minimum event duration is found to be 15 yr, being twice of the minimum 
duration for the distant source case. With increasing the Galactic longitude /, the expected 
event duration increases toward I ~ 90°, where the tangential velocity is dominated by 
the disk velocity dispersion. Toward the outer Galaxy (I = 90° to 180°), the duration is 
between 55 to 70 yr depending on /, with another local maximum at I = 180°. Note that for 
the disk star case the tangential velocity becomes equal to the disk velocity dispersion at 
both I = 0° and 180°, but the expected event duration (t as t) is slightly different from each 
other, since the mean distance to the lens, and hence the mean angular size of the lens are 
different at I — 0° and 180° (this difference of mean lens distance comes from the density 
distribution in the disk). In any case, the event duration cannot be less than 10 yr for the 
disk source case. Hence the event duration seems too long for practical observation when 
one considers that the lifetime of space astro metric missions is 3 ~ 5 yr at longest. 

Note that this long duration for disk-disk event mainly comes from the assumption 
that the source is no further than the galactic center (the distance is assumed to be 8 kpc 
in this paper). In such a disk-disk event, observer, lens and sources move toward nearly 
same direction, and thus the lens-source proper motion is small, being dominated by the 
disk velocity dispersion. On the other hand, for a source star behind the galactic center, 
the relative proper motion of lens is dominated by the galactic rotation rather than the 
disk velocity dispersion (see figure 2). Thus, if one can observe a distant disk star located 
behind the galactic center, one would be able to detect disk-disk events with relatively short 
duration. However, in practice, this is extremely difficult for astrometric space missions like 
SIM and GAIA because of heavy dust extinction in the galactic plane. To detect such an 
event, an astrometric mission at infrared will be necessary. 

For the bulge-source cases (bulge-bulge and bulge-disk lensing), figure 4 shows that the 
event duration can be practically short. For instance, at the galactic center the expected 
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event duration is 8.1 yr for bulge-disk lensing case, and 2.6 yr for bulge-bulge lensing. 
This short event duration comes from large velocity dispersion of bulge stars. In case of 
bulge-bulge lensing, there is another reason for short event duration, that the lens size itself 
is substantially small because of short lens-source distance (since both lens and source are 
in the bulge). Therefore, if one looks for astrometric microlensing at optical/near infrared 
bands, which are the main observational bands of space astrometric missions, it is practical 
to monitor the bulge sources, rather than disk sources. 



5. Effect of Mass Function and Velocity Dispersion 

When an astrometric microlensing event is detected, the lens mass can be estimated 
from the event duration, because the event duration is proportional to the lens mass in 
astrometric microlensing (see, for instance, equation |18|]). As is mentioned earlier, the 
most interesting aspect of equation ( |T8"D is that for distant source case the factor D ds /D s 
approximately becomes unity, and hence the average event duration becomes independent of 
the lens distance. We note, however, that one cannot uniquely determine the lens mass for 
individual events, since there is no way to know the lens-source impact parameters and/or 
tangential velocity for individual events. Also, disk/bulge stellar masses are not unique 
but have a certain range according to the stellar mass function. Therefore, we need some 
statistical treatments when converting the event duration to the mean lens mass. In this 
section, we consider such effects based on more realistic models of the velocity dispersion 
and stellar mass function in the bulge and disk. 

Instead of equation (|J), here we consider more general form of the event rate given by 

r as t = J J J J Anax^E — v± f(v±) <f>(m) cos S dSdv ± dmdD d . (20) 

Here f(v±) expresses the distribution of lens tangential velocity relative to the microlensing 
tube, and (f>(m) denotes the stellar mass function, which are normalized so that 
/ f( v ±) dv± = 1 and / <p(m) dm = 1, respectively. From the relation between v± and £ ast 
(equation |TB|), one can obtain 

dv± = _ 2 P™ xRe cos5<it ast . (21) 
Together with above equations, the event rate may be written as 

r -t = 2 -W~ I I I l P -S- {^) f * cos 3 5 d5 dt as t dm dD d (22) 
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By differentiating equation ( p2|) with respect to t ast , one can also obtain the differential 
event rate c£T ast /<it ast , which shows the distribution of observed event rate against event 
duration. 

Figure 7 shows this (ir ast /(itast against t ast for the cases considered in previous sections 
(QSO-disk, QSO-bulge, disk-disk, bulge-disk, and bulge-bulge lensing). All of the results 
presented in figure 7 are calculated for (I, 6) = (0°, 0°). For QSO-disk and QSO-bulge cases, 
there is a sharp peak at short duration (~ a few years), which corresponds to events caused 
by the stars located at the lower end of the stellar mass function. Thus, if one can trace this 
peaks, it is possible to determine 'typical' lens mass, which is likely to be (or at least fairly 
close to) the minimum cutoff mass of the mass function (provided a sharp cutoff of the 
stellar mass function really exists). For instance, for QSO-disk lensing cases the tangential 
velocity of lens is likely to be around 440 km/s (see figure 2) and the event duration 
becomes independent of the lens distance D& (equation [18|]) . Thus the strong peak of 
the differential event rate is likely to be caused by lenses with rather unique parameters; 
m ~ m-L and v± = 440 km/s. For bulge-disk and bulge-bulge cases, there still exist a peak 
like QSO-disk/bilge lensing cases, but at t ast ~ 0. Also, the probability of lensing is lower 
by more than an order of magnitude. Hence, it is not easy to estimate simply a 'typical' 
lens mass for bulge source cases but one has to take a careful statistical treatment. For 
disk-disk cases, there is no such a sharp peak in <ir ast /(it as t plot. Possible reason for this is 
that for disk-disk case the lens distance widely ranges {{D^/Ds) from to unity) and also 
the tangential velocity varies significantly with lens distance (see figure 2), which causes 
nearly equal-weight contribution of lens with wide range of v± and 

Figure 8, which is a plot of cumulative event rate vs event duration, shows this 
situation from a different point of view. In figure 8, for QSO-disk and bulge-bulge events, 
the cumulative event rate is steeply rising with increasing even duration, demonstrating 
that more than 80% of events have event duration less than 4 yr. For bulge-disk and 
QSO-bulge lensing cases, the rise at short event duration is relatively gentle, but yet 80% of 
events have duration less than 8 yr. In contrast to these cases, the cumulative event rate for 
disk-disk lensing case is only slowly increasing with event duration, indicating considerable 
contributions of events with long duration. 

Note that for QSO-disk and bulge-bulge lensing cases the steep rise in the cumulative 
event rate strongly depends on whether a sharp truncation of the stellar mass function really 
exists. If the stellar mass function has a turnover at a certain mass and gently declining 
toward its lower end, then the rise of the cumulative event rate would be softened. In that 
case, it is not so easy to estimate of a typical mass of 'dominant' population. Alternatively, 
one needs to calculate many of plots like figure 8 with various model parameters as 
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templates, and a comparison between these templates and observed event durations should 
be made to extract reliable information on the lens mass (and possibly the stella mass 
function) from astrometric microlensing events. Thus, figure 8 gives an example of such 
templates to be compared with observed astrometric events. 

6. Discussion and Conclusion 

In this section, we discuss the implication of astrometric microlensing for VERA 
by focusing on the strategy of practical observations. First, here we briefly summarize 
advantages as well as disadvantages of using VERA for astrometric microlensing search. 
Advantages of the distant source case (which is the case for VERA) can be summarized as 
follows: I) large event probability expected from the optical depth and event rate, II) short 
event duration compared to the mission lifetime, and III) less uncertainty in lens mass 
determination, as the expected event duration is independent of the lens distance. On the 
other hand, disadvantages of distant source observation with VERA are: I) small number of 
sources compared to stars in the Galaxy, and II) possible structural variation of the sources 
which could cause an apparent position shift of sources. 

As for the first disadvantage (shortage of sources), it is true that the number of radio 
sources is fairly small compared to the number of visible stars in the Galaxy. For instance, 
the number of compact radio sources that are expected to be observable with VERA is 
around 2000 (e.g., Ma et al. 1998, Peck & Beasley 1998). Moreover, the number density 
of compact sources in the Galactic plane is smaller than that in the off-plane region, since 
surveys of compact sources are conducted mainly in the off-plane regions (e.g., Patnaik et 
al. 1992; Peck & Beasley 1998). Thus, in order to have a sufficient number of sources, one 
has to conduct a survey in the Galactic plane region (e.g., Honma et al. 2000). 

However, even at this stage there exist some sources which can be used for astrometric 
microlensing search. For instance, there is a radio source only 0.7° away from the Galactic 
center (Backer & Sramek 1999; Reid et al. 1999). Also, according to Lazio & Cordes (1998) 
there are at least five extra-galactic radio sources within a few degrees from the Galactic 
center (including the one 0.7° away from the Galactic center). Properties of the five sources 
are listed in table 3, including event rates for disk lens and bulge lens cases. The event rates 
for disk lens cases vary from 6.0 x 10~ 3 to 1.1 x 10~ 2 event/yr, and the total event rate 
for the five sources is found to be 4.0 x 10~ 2 event/yr. Hence, even if one monitors only 
these five sources, one would detect an astrometric microlensing event caused by a disk 
star within 25 yr. The event rate for bulge lens is also quite high, varying from 2.8 to 4.7 
xlO -3 event per year. If the bulge event rate is included, the total event rate for the five 
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sources is 6.0 x 10~ 2 event/yr, indicating that within 16 years at least one of the sources is 
being lensed by bulge or disk star. Note that these five sources are located in about 2° x 8° 
area around the Galactic center (see figure 2 in Lazio & Cordes 1998). If the same number 
density applies to the rest of the Galactic plane, one can expect to find fairly large number 
of radio sources (i.e., a few tens). Therefore, the shortage of sources can be resolved by 
further search for radio sources in the Galactic plane. 

As to the second disadvantage (QSO structure variation), it is not easy to distinguish 
the structure effect and astrometric microlensing within a short period. However, the image 
trajectory of an astrometric microlensing is a perfect circle for most cases except for events 
with extremely small impact parameter. On the other hand, it is unlikely that an image 
motion caused by the structure variation becomes a perfect circle. Hence, if one monitors 
the source position for the full event duration, one can easily discriminate the structural 
effect and astrometric microlensing (but of course there is no way to separate both effects 
if the two effects simultaneously happen to one source, and to avoid this it is better not to 
observe sources with significant structure variation). 

In conclusion, astrometric microlensing events due to stars are practically detectable 
based on the observation of distant radio sources with VERA. If such events are detected, 
one can estimate the lens mass, and this will probably brings us new information on the 
stellar mass function at the lower end. According to the event rate calculation presented in 
the previous sections, it is not easy to detect an event, but it is possible if one monitors a 
few tens of sources for a decade or so. Although this is time-consuming, it is a worthwhile 
study when one considers its scientific importance. 
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Table 1. Summary of the mass model of the Galaxy used in this paper. 



component 


parameter 


value 


bulge 


total mass (M b ) 


8xl0 9 M Q 




scale length (a) 


1 kpc 


disk 


local density (pa) 


O.O8M /pc 3 




scale length (of) 


3.5 kpc 




scale height (h) 


300 pc 




velocity dispersion (era) 


20 km/s 
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Table 2. Summary of source and lens pair considered in this paper. 



name in text 


source 


lens 


A 


QSO-disk 


QSO 


disk star 


oo 


QSO-bulge 


QSO 


bulge star 


oo 


disk-disk 


disk star 


disk star 


8 kpc 


bulge-disk 


bulge star 


disk star 


10 kpc 


bulge-bulge 


bulge star 


bulge star 


10 kpc 
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Table 3. Possible targets for astrometric microlensing search with VERA. 



source name a / b r disk (event /yr) r bulge (event /yr) 



TXS 1741-312 


357. 


.865 


-0.996 


6.0 


X 


10" 


-3 


3.9 


X 


10" 


-3 


TXS 1740-309 


358. 


002 


-0.637 


7.5 


X 


10" 


-3 


4.1 


X 


io- 


-3 


TXS 1739-297 


358. 


.917 


+0.072 


1.1 


X 


10" 


-2 


4.7 


X 


io- 


-3 




358. 


.983 


+0.580 


7.9 


X 


io- 


-3 


4.7 


X 


io- 


-3 


TXS 1748-253 


3.745 


+0.635 


7.3 


X 


10" 


-3 


2.8 


X 


10" 


-3 



total (E T) 4.0 x 10- 2 2.0 x 10" 2 



a Source name in the Texas survey (Douglas et al. 1996). In Lazio & Cordes (1998) 
source names are given by the galactic longitude and latitude. 
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Fig. I. — Optical depth for QSO-disk, QSO-bulge, disk-disk, bulge-disk, and bulge-bulge 
lensing cases. Thin lines are optical depth for QSO-disk lensing for b = 0°, b = 2.5°, and 
b = 5°, respectively. Dashed lines are for disk-disk lensing (for b = 0°, and 6 = 5°). Dotted 
lines are for QSO-bulge lensing, for bulge-disk lensing, and for bulge-bulge lensing cases from 
top to bottom. 
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Fig. 2. — Distribution of tangential velocity of the lens relative to microlensing tube for 
I = 0°, 10°, 20°, and 30° cases (thin lines for QSO-disk lensing cases, and dashed lines for 
disk-disk lensing cases). For QSO-disk lensing, the tangential velocity becomes close to 2v ciTC 
(~ 440 km/s) behind the Galactic center. 
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Fig. 3. — Distribution of event rate with galactic longitude / (notations are the same to 
figure 1). For QSO-disk lensing (thin lines), three cases for b (0°, 2.5°, and 5°) are shown. 
While the event rate for QSO-disk/bulge lensing case becomes maximum toward the Galactic 
center, the event rate for disk-disk lensing (dashed lines) is largest around / = 10°. 
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Fig. 4. — Distribution of expected event duration ((t as t)) for b = 0°, b = 2.5°, and b = 5°, for 
both QSO-disk (thin lines) and disk-disk lensing cases (dashed lines). Also shown are event 
duration distributions for QSO-bulge lensing, bulge-disk lensing, and bulge-bulge lensing 
(dotted lines, shown only for / within 20°). 
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Fig. 5. — Distribution of differential event rate with lens distance for QSO-disk lensing case 
(for / = 0° and 10°). The differential event rate sharply peaked behind the galactic center. 
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Fig. 6. — Distribution of averaged event duration with lens distance for QSO-disk lensing 
case (for I = 0° and 10°). The disk stars behind the galactic center cause events with short 
duration (~ 4 yr), while those near the Sun cause much longer event (~ 90 yr). 
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Fig. 7. — Distribution of differential event rate (c£r ast /t a st) against event duration t ast (for 
I = 0° and 6 = 0°). For QSO-disk/bulge lensing cases, the differential event rate has a sharp 
peak at short duration of a few years, while that for disk-disk lensing is rather flat. 
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Fig. 8. — Cumulative event rate against event duration for QSO-disk, QSO-bulge, disk-disk, 
bulge-bulge and bulge-disk lensing cases. QSO-disk and bulge-bulge cases show steep rises 
of cumulative event rate at short event duration, in contrast to disk-disk lensing case. 



